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Abstract 

A path integral evaluation of the Green's function for the hydrogen atom ini- 
tiated by Duru and Kleinert is studied by recognizing it as a special case of the 
general treatment of the separable Hamiltonian of Liouville-type. The basic dy- 
namical principle involved is identified as the Jacobi's principle of least action for 
given energy which is reparametrization invariant, and thus the appearance of a 
gauge freedom is naturally understood. The separation of variables in operator 
formalism corresponds to a choice of gauge in path integral, and the Green's func- 
tion is shown to be gauge independent if the operator ordering is properly taken 
into account. Unlike the conventional Feynman path integral, which deals with a 
space-time picture of particle motion, the path integral on the basis of the Jacobi's 
principle sums over orbits in space. We illustrate these properties by evaluating 
an exact path integral of the Green's function for the hydrogen atom in parabolic 
coordinates, and thus avoiding the use of the Kustaanheimo-Stiefel transformation. 
In the present formulation , the Hamiltonian for Stark effect is converted to the one 
for anharmonic oscillators with an unstable quartic coupling. We also study the 
hydrogen atom path integral from a view point of one-dimensional quantum grav- 
ity coupled to matter fields representing the electron coordinates. A simple BRST 
analysis of the problem with an evaluation of Weyl anomaly is presented . 
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1 Introduction 



The Hamiltonian of the hydrogen atom provides one of those few examples which can be 
solved exactly, and as such various interesting alternative methods to solve it have been 
proposed in the past. Among those methods, one may count an elegant momentum space 
analysis by Schwinger[l]who exploited the 0(4) symmetry hidden in the hydrogen atom 
Hamiltonian. The path integral treatment of the problem is also interesting not only for 
a methodological interest but also for a pedagogical purpose. A comprehensive study of 
a semi-classical approximation of the path integral for the Green's function at a given 
energy has been performed by Gutzwiller[2]; in particular, he found an exact Green's 
function for negative energy in the polar coordinates of momentum space by a semi- 
classical approximation. It is also known that a path integral of the s-wave propagator ( 
or evolution operator ) is obtained by summing a perturbative series [3]. In 1979, Duru 
and Kleinert[4] showed an elegant path integral method to evaluate the Green's function 
for the hydrogen atom exactly. Two basic ingredients in their method are the use of 
a re-scaled time variable and the so called Kustaanheimo-Stiefel transformation[5]which 
reveals the 0(4) symmetry explicitly in the coordinate space. Many of the clarifying 
works of this approach have been published [6] - [15]. The main issue in these works 
is a physical meaning of the "re-scaled time variable". We here study this issue from a 
completely different view point on the basis of the Jacobi's principle of least action by 
using a general gauge theoretical technique. 

We first recognize the procedure in Ref.[4] as a special case of the general treatment 
of classically separable Hamiltonian of Liouville-type. The basic dynamical principle in- 
volved is then identified as the Jacobi's principle of least action for given energy. The 
fundamental feature of the conventional Feynman path integral, which is associated with 
the Hamilton's principle of stationary action, is that it deals with a space-time picture of 
particle motion. On the other hand, the path integral on the basis of the Jacobi's prin- 
ciple of least action is basically static and analogous to geometrical optics. A space-time 
picture is thus lost, and one deals with a sum over orbits in space instead of space-time. 
Another characteristics of the Jacobi's principle of least action is that it is reparametriza- 
tion invariant, and the appearance of a gauge freedom to fix an arbitrary parameter, 
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which describes orbits for fixed energy, is clearly seen. The general technique of gauge 
theory is thus applicable to the evaluation of path integral, and a suitable choice of gauge 
simplifies the problem such as the hydrogen atom. In particular, the Green's function is 
shown to be gauge independent. 

The use of the Kustaanheimo-Stiefel transformation is rather technical and it is not 
essential in solving the problem exactly. We in fact show a simple trick in parabolic 
coordinates which solves the hydrogen atom exactly. This trick was used before in a 
different context by Ravndal and Toyoda[16]. We also note that the path integral in 
holomorphic variables (or the coherent-state path integral) is convenient to evaluate the 
path integral of the hydrogen atom. 

As an application of the present approach, the Hamiltonian of Stark effect is shown 
to be reduced to anharmonic oscillators with an unstable quartic coupling, for which a 
resummation technique of a perturbative series is well-known. We also study the hydrogen 
atom path integral from a view point of one-dimensional quantum gravity, and a simple 
BRST analysis of the problem with an evaluation of the Weyl anomaly is presented. This 
provides an alternative way to see the gauge independence of Green's functions. 

2 Separable Hamiltonian of Liouville-type 

We here explain the basic procedure to treat a general separable Hamiltonian of Liouville- 
type. We start with a separable Hamiltonian 

H = v( * , A ^-(p\ + p\) + UM + U 2 (q 2 )} (2.1) 

where the variables change over 00 > q±,q 2 > —00. A general Hamiltonian of Liouville- 
type is given by 

H = v^Uw^ ^m^ + + Ul(Ql) + u * m) (2 - 2) 

but after a canonical transformation 

1 „ rQ 



Pi=Pi , / JW 1 (Q)dQ = q 1 



1 1 P2=P2 , \/W 2 (Q)dQ = q 2 (2.3) 



and a suitable redefinition of V and U, we can derive a Hamiltonian of the form in (2.1). 
We may then solve the Schroedinger problem 

Ei> = vmIvm ^*® + U > M + u ^ (2 - 4) 

where 

p t = -ih-£- (2.5) 
dqi 

for 1 = 1, 2, and the volume element dV, which renders the Hamiltonian H in (2.4) hermi- 
tian, is given by 

dV=(V 1 (q 1 ) + V 2 (q 2 ))dq 1 dq 2 (2.6) 

The classical Hamiltonian (2.1) does not completely specify the operator ordering in (2.4), 
and the simplest ordering is adopted here. A precise operator ordering needs to be fixed 
depending on each explicit example; a concrete example shall be given for the hydrogen 
atom later. 

One may rewrite the above Schroedinger equation (2.4) as 

H T ip = (2.7) 
with a total Hamiltonian defined by a specific gauge condition, 

Ht = ^(Pl+Pl) + U M) + U 2 {q2) ~ E{y-y{gx) + V 2 {q 2 )) (2.8) 

The meaning of a total Hamiltonian is clarified later. A general procedure to deal with a 
completely separated operator Ht is to consider an evolution operator for a parameter r 
defined by 

(9i6, q2 b \e- lflTT/h \q la , q 2a ) = (qib\ exp [-(i/h)(^-pj + U^q,) - EV 1 (q 1 ))r}\q la } 

x (q 2b \ exp [_(i//i)(J-p2 + U 2 (q 2 ) - EV 2 (q 2 ))r]\q 2a ) 
= J pg 1 pp ie ( l / n )/o T fr l9 ' 1 -(2kP?+ l/l (« 1 )- £yi ( < ? 1 ))> dr 
x y > pg 2 pp 2 e (i/a) /o T ^ 2 9' 2 -(^ +% ( ( ' 2 )- W2 ( <i,2 ))> dr (2.9) 

The parameter r is arbitrary , and by integrating over r from to oo one obtains a 
physically meaningful quantity 

(qib,q2b\-~—\qi a ,q2a) semi— classical 



H 



T 



i / dr- 
Jo 



pKih{d qi {T)/d Pi m qia ^2mh(dq 2 (r)/dp 2 (0)) q2a 
x exp {(i/h)S cl (q lb , q la , r) + (i/h)S c fa b , q 2a , r)} (2.10) 

where we wrote the result of a semi-classical approximation for the path integral [17] [18] 
[19], though in certain cases one may be able to perform an exact path integral in (2.9). 
The pre-factor in (2.10) is written in terms of classical paths , for example, 

qici(r)=q l (r- 1 q la ,p 1 (0)) (2.11) 

Namely, the classical paths dictated by the total Hamiltonian H T are expressed as func- 
tions of the initial positions and momenta. On the other hand, the classical action S d 
is expressed as a function of the initial position, final position and elapsed "time" r by 
eliminating pi(0) dependence;for example, 

S d (q u , qia, r) = F { Pl q\ - {^-p\ + Ufa) - EVfa))} d dr (2.12) 
Jo 2m 

with ?i(t) = qib- If one solves the Hamilton- Jacobi equation in the form 

S{qu, Qia, t) = -At + S(q lb , q la ; A) (2.13) 

one treats A as a dynamical variable and regards the above equation as a Legendre 
transformation defined by 

dS(q lb ,q la ;r) 



dr 

dS(q lb ,q la ;A) 



= -A 

= t (2.14) 



OA 

The variable A is then eliminated. This may be regarded as a classical analogue of 
uncertainty relation; if one specifies r, the conjugate variable A becomes implicit. It is 
known that the semi-classical approximation (2.10) is exact for a quadratic system such 
as a simple harmonic oscillator [18]. 

We next note the relation for the quantity defined in the left-hand side of (2.9) 

(?16,?2b|-£-kla,?2a) 

I 1 1 

— \qib,q2b\- j \Qla, <72a/T77 \ . \r ( \ 

W)W gr vM + vfa a ) 

= {qib,q2b\-^ -\qia,q2a) v , r—TTJ- s ( 2 - 15 ) 

H-E Vi{qi a ) + V 2 {q 2a ) 

= 777 h~~B \ F?i / — 7~ ; ; — ~ {qib, q2b\qia, q2a) i ==} 

H (^,^ b ,..)-E ^Vfa b ) + V 2 {q 2b ) y/vfal) + Vfaa) 



by recalling (AB) 1 = B l A 1 . The state vectors in these relations are defined for the 
volume element dq\dq 2 as 

J dq 1 dq 2 \q 1 , q 2 ) (q 1 , q 2 \ = 1 (2.16) 

Note that the definition of the 5-function in (q[, q' 2 \qi, q 2 ) = S(q[ — qi)5(q' 2 — q 2 ) depends 
on the choice of the volume elememt in (2.16) and thus on the choice of Ht- The last 
expression in (2.15) is thus correctly defined for the original Hamiltonian H and the orig- 
inal state ip in (2.4) with the volume element dV in (2.6), since we have the completeness 
relation from (2.16) 

f dV 
J Vi(gi) + V 2 (q 2 ) 

The left-hand side of (2.15) thus defines the correct Green's function for the original 

operator (H — E)~ l by noting the symmetry in q a and q b . One can then define the 

conventional evolution operator by 

(gib, q2b\e- iH{tb - ta)/h \q la , q2a)conv 

I POO I 

= ^-r\ dEe-^-^\q lb ,q 2b \ J <h;<ha) v( r~ „ , , (2.18) 

2mhJ-oo H_ le -E Vi(qi a ) + V 2 (q 2a ) 

where e is an infinitesimal positive number. The total Hamiltonian changes for a different 

choice of gauge condition in the Jacobi's principle of least action to be explained below. 

Consequently, the volume element, which renders H T hermitian, generally depends on the 

choice of gauge. In this case, one has the relation 

1 



(qib, q2b\qia, q2a) conv ~ (?16, <?2fc | <?la, <?2a) 



Vi(q la ) + V 2 (q 2a ) 

--(qib,q2b\qia,q2a) 



'K(g 16 ) + y 2 ( g26 ) " ^(q^ + V^) 

The meaning of the total Hamiltonian Ht (2.8) becomes transparent if one starts with 
the Jacobi's principle of least action for a given E 

S = f drL = f T dr^miEiV^) + V 2 (q 2 )) - (C^fo) + U 2 (q 2 ))}(q! + (£) 
Jo Jo v 

= / ^2m[E(V 1 (q 1 ) + V 2 (q 2 )) - (U 1 (q 1 ) + U 2 {q 2 ))}[{d qi f + {dq 2 f} (2.19) 

which is reparametrization invariant. One then defines the momenta conjugate to coor- 
dinates 

Pi = ttt = v/2m[E(K(gi) + V 2 {q 2 )) - (U 1 (q 1 ) + U 2 {q 2 ))} x J (2.20) 
°Ql yj(q{ + g 2 2 ) 



and obtains a vanishing Hamiltonian, which is a result of reparametrization invariance, 
and a first class constraint as the generator of reparametrization gauge symmetry, 

H = p,g,-L = 

<f>{<luPi) = v( * , A ^-ipl+pD + UM + UMj-E^O (2.21) 
+ V 2 {q2) 2m 

Following Dirac[20], one may then define a total Hamiltonian 

H T = H + a(q h pi)(j)(q h pi) 

= a(5i,Pi)0(«,Pi)-O ( 2 -22) 

where an arbitrary function a(qi,pi) specifies a choice of gauge or a choice of the arbitrary 
parameter r in (2.19), which parametrizes the orbit for a given E. The quantum theory 
is defined by( up to an operator ordering) 

d 

ih—^j = H T ij (2.23) 
or 

with a physical state condition 

&(oi,Pi)4>(gi,Pi)^phv = Q (2-24) 

A choice of the specific gauge ct(qi,pi) = V\(qi) + V^t^) gives rise to (2.7) and the choice 
a (<lhPi) — 1 gives the conventional static Schroedinger equation (2.4), since ip appearing 
in these equations are physical states. 

The basic dynamical principle involved is thus identified as the Jacobi's principle of 
least action, which is analogous to geometrical optics, and the formula of an evolution 
operator (2.9) dictated by (2.23) provides a basis for the path integral approach to a 
general separable Hamiltonian of Liouville-type. The path integral in (2.9) deals with 
a sum over orbits in space instead of space-time, and the notion of re-scaled time does 
not explicitly appear in the present approach; the evolution operator (2.9) essentially 
generates a gauge transformation. 

3 Hydrogen Atom 

3.1, Analysis in Parabolic Coordinates 
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We analyze the hydrogen atom by starting with the Hamiltonian written in terms of 
parabolic coordinates 



where the parabolic coordinates are introduced via the cylindrical coordinates (p, ip, z) by 



v = ^y + z) 



VP 2 + ^ (3-2) 



and ip stands for the azimuthal angle around the z axis. We further perform a canonical 
transformation which simplifies the kinetic term in H as 

1 2 

C — ~U , < W < OO 

4 

>J&t = Vu 

1 2 

rt — -v , < t> < oo 

4 

v 7 ^ = Pv (3.3) 
and the Hamiltonian becomes 

H = ^-(-A—i) \vl + — 2 V% + Vl + \v\\ - -A-^ e2 (3-4) 
2m u z + v z u v z * u z + v z 

where r = £ + r/ = (u 2 + t> 2 )/4. This Hamiltonian is not yet a separable one of Liouville- 
type. 

One may solve the Schroedinger equation 

Hip = Eip (3.5) 

or equivalently 

H T i) = (3.6) 

with 



where uj is defined by 



I mtJ 2 = --E (3.8) 
2 4 v ; 
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We consider the case E < for the moment. Ht stands for the total Hamiltonian defined 
by a specific gauge condition; a general definition of Ht will be given later in (3.47). 
Eq.(3.6) may be rewritten in an equivalent form as 

H T ij = 

(p v -Pv>)i> = (3.9) 
We here introduced auxiliary variables (p^',^') as 

Ht = ^iPl + ^+pl + ^ + ^i^ + v^-e 2 

and we defined 

u = (ui, u 2 ) = (u cos ip, u sin ip) 

Ub 



v = {y i, V2) = (v cos <f', v sin ip') 

$ = Pl + ^Pl' (3-11) 



IT 

The subsidiary condition in (3.9) replaces the use of the Kustaanheimo-Stiefel transforma- 
tion, and at the same time it renders a Hamiltonian of Liouville-type. This introduction 
of auxiliary variables (3.11) has been discussed by Ravndal and Toyoda[16]. The use of 
the subsidiary condition (3.9) in place of the Kustaanheimo-Stiefel transformation may 
be useful for a pedagogical purpose. 

A general procedure to deal with a completely separated operator such as H T in (3.10) 
is to consider an evolution operator for a parameter r defined by 

(u b ,v b \e-^ n \u a ,v a ) 

= e^(u b \exp[-(i/h)(^-f u + ^fu 2 )r]\u a ) 
x(v b \exp[-(i/h)(^f v + ^-v 2 )r}\v a ) 

= c »e 2 r/ mUJ )4/2 

2irih sin ujt 

xexp{—-. \(u 2 b + vl + u 2 a + v 2 a ) cosujt - 2u b u a - 2v b v a ]} (3.12) 
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where we used the exact result for a simple harmonic oscillator [18] 

(q b \exp[-{i/h)(^-p 2 + ^-q 2 )T}\q a ) 

, mw .1/0 , imuj r , 9 9 . n .„ „ „. 

= — ) ' ex p{^r~. — M + ?« cos ^ - 2 wa } 3.13 

27rmsmcjr 2ft sm cur 

which can be established either in path integral or in operator formalism( see (2.10)). 
Also, the harmonic oscillator is symmetric with respect to coordinate or momentum space 
representation [2]. 

A crucial observation here is that p v and p^ are preserved during the evolution dictated 
by the operator Ht in (3.10), since Ht] = [jy, Ht] — 0. It is then sufficient to impose 
the constraint (3.9) only on the initial state , for example. Starting with a general state 
belonging to the eigenvalues p v = m and p^ = m! 



e , " , V mV (3.14) 

we can use the following trick 

J** ^. e im( v +e) e im'( v '-e) = 6m ^ e i m ^+v') ( 3 . 15 ) 

to project out the state satisfying p v = p v r, and y? + y?' is regarded as the actual azimuthal 
angle. 

We thus obtain 

(u b , v b , (<p + V ') b \e- iHTT/h \u a , v a , (ip + <p') a ) 
iJi^, x9 f 27T d6 r imui 



/ — exp{ -— [(u b + v b +u a + v a ) cos ujt - 2u b u a - 2v b v a \ \ 

Jo In sin cur 



2mh sin cur Jo 2n 2h sin cur 
ie 2 T mu 2 

= e te T ( y x 

2-Kin sin ujt 

r 2n d6 imuo r . _ ^ x1 / 9 

/ 7r ex P\7^-- 4 (&* + ?6 + Va + Vb) cos cur - AV2(r a r b + x a x b ) 1 cos(6> + 7) } 

Jo 2n 2n sm ujt 

ie 2 T , muj 2 r 2imcu . . 2\plimuj _ _ 1/2 

= e ( o •* • ) ea; P{"r^ ( r « + rfe ) coscur}/ (-r- f fl r t + x a x b ) 1 ) 

2-nihsmujT asm ujt asm ujt 

(3.16) 

In this evaluation we start with the relation 

u b u a + v b v a = u b u a cos Aip + v b v a cos Ay?' (3-17) 



10 



with A(p = <p b - <p a , Atp' = (p' b - (p' a , and 



u b u a cos(A(/9 + 9) + v b v a cos(A<// — 9) 
= {u b u a cos Aip + v b v a cos A<//) cos 9 
+ (—u b u a sin Ay? + v b v a sin Ay/) sin 9 



4\/6>£a + Wa + tbtaVbVa cos(Ay? + Ay?') cos(6» + 7) 



= 2\/2 y r a r b + z a z b + p a p 6 cos(Ay? + A</?') cos(# + 7) 

= 2y/2y/r a r b + x a x b cos(# + 7) (3.18) 

where we used the definition of variables in (3.3), and 7 is a number independent of 9. 
We also defined a modified Bessel function 

2V2imu _ 1/2 /"^ d9 2^/2irmj(r a r b + f a ^b) 1/2 ^ ,„ in x 
7 ^ r a r 6 + a; a a; b 7 = / — exp{ — cos#} 3.19 

The parameter r is arbitrary, and we eliminate r to obtain a physically meaningful 
quantity by 



i / dr(u b , v b , (cp + ip') b \e~ lHTT/ \u a , v a , (ip + <p') a ) 
Jo 

= (u b ,V b , (if + ip')b\-*r-\u a ,V a , (ip + <p') a ) 

. f°° , J> 2_, muo 2 2?mw , . , 

= 1 dre ( , ) exp{— (r a + r 6 ) coswr} 

Jo zm/isiiiwr nsmujT 

T ,2y/2imw , _ _ s1 /2 . 

ftsmcjr 



ma; r 1 ,-2muu . I -\- A Ay Imu A' ^^ xl / 2x 

« Jo dXX TT^W exp[ — {ra + n)i —\ )]Ioi ^^—\ {ran + XaXb) ] 



(3.20) 



where we rotated r by 90 degrees , r — > —it, and defined 

A = e~ 2u)T 

v = e 2 /2co (3.21) 

We next show that (3.20) gives an exact Green's function for the hydrogen atom by 
noting the sequence 

(u b ,V b ,(p b \^-\u a ,V a , (f a ) 
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(tb,Vb,(Pb\ -, I -\U,Va,(fa)( 



= (Zb,r]b,<Pb\ £■ h -\£a,ria,<Pa)( -Z ]_ ) 

H(£, T],ip) — E £a + Va 

ai^aifa) rp ■ ) 



f T// . \ n \ /t= — ; \S0> '/o> ^o sm '/a, rr. — ■ , 

H(£ b , Vb, <Pb) - E V?6 + V £a + ??a 

= i<^b l^a) (3.22) 

47r p /2m - e 2 /r - £ 

where we used f in place of f + <// and the relation (AB) -1 = B^Ar 1 . 
The volume element changes in this transition from to i7 as 

dV = 2nuvdudvdf 
-^dV = (f + rj)dV = 4tt x 2(f + rj)d^drjdf 

= 4tt x r 2 drd cos Odip (3.23) 

The bra- and ket- vectors in (3.22) are normalized in the combination 

J dVo\u, v , <f) (u, v , ip\ = 1 
/ dV\^,r],if)-^—(^,r],if\ = 1 

J d 3 x\x){x\ = 1 (3.24) 

and the extra factor of An in dV = Anr 2 drd cos ddip requires the appearance of the factor 
of 1/47T in the last expression in (3.22). The appearance of 2ir in dV is an artifact of the 
variable if' in (3.11). This normalization condition of bra- and ket- vectors together with 
a symmetry in x a and Xb justify the identification (3.22). A more explicit and concrete 
analysis of eqs.(3.22)~ (3.24) will be given in connection with the Jacobi's principle later. 
As for the operator ordering, the momentum operator changes in (3.22) as 

fu+Pl = (-) 2 [-d u ud u + -d v vd v ] 

% U V 

= Pz&t + PnVPr, (3-25) 



and 



1 1-2 

)(p&Pt + PriVPv) + J^Pl = P (3-26) 
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where the right-hand side is written in cartesian coordinates. We note that dV® and dV 
in (3.23) respectively render and H(£,r),(p) hermitian. 

Combining (3. 20), (3. 22) and (3.26), we have thus established the exact Green's func- 
tion including the operator ordering 

._. h ._. 2m 2 u r 1 1 -2mu . W l + A,. 

{x b \— - \x a ) = —2- / dXX - —exp[— — (r a + r b ){- -) 

p /2m - e 2 /r - E nh Jo (1 - A) 2 h 1 - A 

. - A\f2mw X 1/2 , _ ^ Nl/9N /„ 

x/ (^ — Y^x^ aTb + a h) ] ( 3 7) 

It is known that this formula, which was first derived by Duru and Kleinert[4], is a 
Fourier transform of Schwinger's momentum space representation [1]. The continuation 
to the scattering problem with E > is performed by the replacement 

iq — > v — > iv (3.28) 

in the above formula. 

One can understand the spectrum of the hydrogen atom by looking at H T in (3.10). 
This problem, which has been analyzed by Ravndal and Toyoda[16], is briefly summarized 
here in connection with the Jacobi's principle of least action and the Stark effect to be 
discussed later. If one defines the oscillator variables 



a k = - 7 =[\j—u k + = Pu k \, 



one obtains 



1 Imu 

~a k = _[^_ t , fc + _ A , fc ], A; = 1,2 (3.29) 



2 

,t„. , _L 01 _ .2 



H T = hu[J2(% a k + %a k ) + 2] - e 2 
fe=i 

p v = ih[a\a 2 — c4 a i] 

Ptp' = ih[d\(i2 — a\ai\ (3.30) 



After a unitary transformation 



di = -^(61-162) 

a 2 = ■^=(-ib 1 + b 2 ) (3.31) 
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and a similar transformation of a\ and a 2 , one obtains 

2 



H T = hcu[J2(b{b k + b{b k )+2]-e 2 
k=l 

p v = hlblh-blh] 

Pip , = hiblk-blk] (3.32) 



By defining the number operators 



nk = b\b k , 

n k = blh, k = 1,2 (3.33) 

the total Hamiltonian is given by 

H T = hcu[ni + n 2 + hi + h 2 + 2] — e 2 

= %uo[2ni —p v /% + 2n 2 + p (p >/h + 2] - e 2 

= 2^[nx + n 2 + 1] - e 2 (3.34) 

by noting the physical state condition p v = p^i. 

We thus define the principal quantum number (or its operator) n by 

n = m +h 2 + 1 = 1,2,3, (3.35) 

and the physical state condition 

(2nhuj — e 2 )ip p hy S = (3.36) 

gives rise to the Bohr spectrum 

by noting the definition of u in (3.8). 

As for the degeneracy of states with a fixed n, we have n combinations of (ni,n 2 ) 
given by 

n 2 = n — (ni + 1), 

ni = 0,1, ,n-l (3.38) 
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For each fixed (ni,n 2 ), we have the constraints 



^2 + Pip'/h — n — (ni + l)+ p^'/h > 



(3.39) 



which gives 



n\ > p^/h = p^'/h > —n + (ni + 1) 



(3.40) 



Namely, we have 



ni — (—n + ni + 1) + 1 = n 



(3.41) 



possible values of p v . We thus have 



n x n = n 2 



(3.42) 



degeneracy of states with a fixed principal quantum number n, as required. This somewhat 
unorthodox classification of states may be useful in the analysis of the Stark effect to be 
commented on later. 

Incidentally, the formula (3.32) suggests that one can define the Green's function by 



for a direct path integral in (3.16), instead of using the trick noted in (3.15). The 9 integral 
imposes the constraint p^ = p^. The path integral of an evolution operator in terms of 
oscillator variables is known as holomorphic ( or coherent state) path integral[21]. It is 
interesting that (3.43) resembles a closed string propagation on a cylindrical world-sheet. 
3.2, Jacobi's Principle of Least Action 

The meaning of the total Hamiltonian in (3.7) becomes transparent , if one starts 
with a Nambu- Goto- type Lagrangian ( the Jacobi's principle of least action for a given E 
) which is reparametrization invariant, 




fc=i 



(3.43) 



S 




(3.44) 



The momenta conjugate to coordinates are then defined by 




(3.45) 



15 



and one obtains a vanishing Hamiltonian as a result of reparametrization invariance and 
a first-class constraint <fi, which is the generator of reparametrization gauge symmetry, 



H = px — L — 

-2 

<P(x,p) = + V(r) - E ~ (3.46) 

Zi ill 

Following Dirac[20], the total Hamiltonian is defined by 

Ht = H + a(x,p)(f)(x,p) 

= a(x,p)<f>(x,p)~0 (3.47) 

and the function a(x,p) specifies a choice of gauge and fixes the arbitrary parameter r 
in (3.44), which parametrizes the orbit for a given E. A change of the parameter r to 
t — 5(3(t,x,p) is generated by 5(3(r,x,p)a(x,p)(f)(x,p) = 5j3(r,x,p)HT, for example, 

Sx(t) = x'(t) — x(t) 

= x(t + 5(3) - x(t) 

= {x,5(3H T } PB 

= 8f3(r,x,p)-^x(r) (3.48) 

in terms of the Poisson bracket, since x'(t — 5(3) = x(t). 
Quantization is performed by 

ih^-i> = H t i/j (3.49) 



with a physical state condition 



dr 



a(x,p)4>(x,p)ij phy = (3.50) 



A specific choice of the gauge a(x,p) = r = £ + 77 leads to the Hamiltonian Ht in (3.7) and 
the choice a(x,p) = 1 gives the original static Schroedinger equation(3.5), since the states 
ip in (3.7) and (3.5) are physical states. Eq.(3.49) gives rise to the evolution operator in 
(3.12), but the parameter r need not be interpreted as a re-scaled time. In fact, the 
evolution operator (3.12), which essentially generates a gauge transformation, deals with 
a sum over orbits in space instead of space-time. 
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We now explain the relations (3.22)~(3.24) in a more concrete manner. We start with 
eq.(3.22) for a generic negative E 

G(E;£b,Vb,<Pb;£a,Va,iPa) = (€b,Vb,<Pb\ ~ — ^ -\^a,Va,<Pa) 

= ^ 4> n {E; rj b , (p b ) 0*(ff ; £ a , rj a , ip a ) (3.51) 

with 

H T (E;^r],ip)(f) n (E;^r],ip) = X n (E)(f) n (E;^,r],ip) 

J <i>n(E;Z,v,<p)ME;Z,v^)dVo = S n:l 

dV = 4tt x 2d^dr]dip (3.52) 

where we used the result in (3.34) and also the variables (£,77,^2) instead of (u, v,<p) 
for notational simplicity. The summation over n in (3.51) is formal including the n 2 
degeneracy. Note that the complete orthonormal states {4> n } in (3.52) are all unphysical 
off- shell states. In path integral, the summation in (3.51) is exactly evaluated in (3.20). 

We next rewrite G(E; £&, rjb, <f b ] (,a, Va, Va) m terms of physical on-shell states by writing 
an unsubtracted dispersion relation (i.e., paying attention only to the pole structure in 
E) as 

G(E; ?7 6 , if b ; 7] a , if a ) = ^ 4>n{E n ] Vb, <Pb) — —7 dx (E U 0„(K; fa, (f a ) 



(E n -E)(- 9 -^) 



(3.53) 



by noting 



X n (E) = \ n {E n ) + (E — E n )^^ = {E- E n ){^r) (3.54) 



for E 

When one defines 



one can show the orthonormality relations 



Jr n (En,ti,V,<p)MEi;ti,V:<p)(ti + v)dVo = Jr n (E n ;^ V ,if)MEr,^V^)dV 

= 5 nil (3.56) 
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with dV = (£ + r))dVo. First of all, from the physical state condition 

H T (E n ;£,r],ip)(f) n (E n ;t,r},ip) (3.57) 

= {^(Pdh + PvVPri) + g^(^ + ^)Pl ~ E n (£ + rj) - e 2 }<p n (E n ] £, 77, v?) = 

one can establish the orthogonality relation 

(E n - Ei) J fa(E n -£,r 1 ,ip) ( f )l (E l ;Z,r 1 ,ip)(Z + r 1 )dV o = (3.58) 

for n ^ I. Also from the relation (3.57) and the fact that the "eigenvalue" e 2 is equally 
distributed for the kinetic and potential terms for harmonic oscillators(in terms of u and 
v variables) , we have 

-E n j ^(E n ,^ V ,ip)MEn^,V^)^ + v)dV = j (3.59) 

namely 

dx n ( En) I ti(K;Z,v,v)ME n ;S,v,<p)(S + v)dVo = i (3.60) 

dE n 

by noting (3.54). This proves (3.56). 



From (3.53), we finally arrive at the expression 

E„ — E 



G(E;£b,Vb,<Pb;ta,Va,<Pa) = ^n(K', €b, Vb, <Pb) F ~ ^ n( E n, £«, Va, <Pa) 



= (tb,Vb,<Pb\-f— -\£a,Va,<Pa) 

= J-(^Ito fia) (3-61) 

47r p I (2m) — e 2 /r — E 

which establishes the gauge independence of the Green's function for negative E. See 
also Ref.(14). Although we here used the same notation for the state |£, 77, <£>) in (3.51) 
and (3.61), the meaning of these states are quite different. This difference is explicitly 
exhibited in (3.22) ~ (3.24). It is important to realize that the exact path integral is 
performed for the off-shell states in (3.51). As for the case of positive energy, one can 
define the Green's function by analytic continuation in the starting expression in (3.51) 
and in the final expression in (3.61). 

More generally, one can establish the gauge independence of the Green's function for 
an arbitrary choice of gauge condition a(£, 77, (p) = /(£, 77, (p) by noting 

H T (E; £, 77, p)(f) n (E; £, 77, <p) = [/(£, 77, <p)H(£, 77, <p) - /(£, 77, ip)E}(f) n (E; £, 77, <p) 

= \ n (E)ME;Z,ri,v) (3-62) 
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Namely 

J <j>* n {E- £, 77, v?)^t(£; £, ?7, v^)0n(^; £, v, vW f = \ n {E) J <p* n (E- £, 77, £, ?7, </?W 

= A n (£) (3.63) 



dV f = — dV (3.64) 



with 

1 

which renders H T (E;£,rj,(p) defined by the gauge a = f(£,rj,ip) hermitian. From (3.62) 
and (3.63), one derives 

= -J <j)* n (E; £, V , V)ME; £, ?7, (3.65) 

If one uses this relation for E = E n in (3.53), one arrive at the expression (3.61) starting 
with an arbitrary gauge condition a = f(£,r),ip) . 

4 A related topic : Stark effect 

As an interesting implication of the present treatment of the separable Hamiltonian of 

Liouville-type, we comment on the Stark effect (the hydrogen atom inside a constant 

external electric field S) described by a Hamiltonian 

1 e 2 

H = p 2 e£z 

2m r 

We thus analyze the total Hamiltonian defined by 

fr 1 -B , mUj2 -2 . 1 -a . ^^2 2 1 /-4 -4n /„,y\ 

^=2^P t , + — « + + — « " e ~- A 9(v -u) (4.2) 

with a constraint 

P v -p' v = (4.3) 
for the coodinates defined in (3.10). The quartic coupling constant g is given by 

g = (4.4) 
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We can thus analyze the Stark effect on the basis of 

(u b ,v b \e- iATT / n \u a ,v a ) 
= e ie "^u b \e-^^ + ^ +1 ^\u a ) 

= e^ j Vp u VueMm) [[Pun - {^f u + ^fu 2 + \g^)]dr} 

x / Vp v Vve W {(m /V^- (±f v + ^fv 2 - \g^)\dr} (4.5) 

This problem in the conventional formulation has been recently analyzed by K. Hi- 
raizumi, Y. Ohshima and H. Suzuki [22] as an application of the resummation technique 
of perturbation series for quantum tunneling [23], which was established for a system of 
anharmonic oscillators [24]. It is interesting that this problem is in fact identical to the 
(two-dimensional) anharmonic oscillator with an unstable quartic coupling in the present 
formulation. 



5 One-dimensional quantum gravity 

A way alternative to (3.44) to see the physical meaning of the parameter r is to study 
the one-dimensional quantum gravity coupled to matter variables x defined by 

[ exp{i f Lhdr} (5.1) 

J gauge volume Jo 

with 

L -^f/- v ^ + E < 5 - 2 > 

where h stands for the einbein, a one-dimensional analogue of vierbein h®, and h = ^fg 
in one-dimension. In this Section, we set h — 1. If one uses the solution of the equation 
of motion for h defined by the Lagrangian C = Lh, the action in (5.1) is reduced to the 
one appearing in the Jacobi's principle of least action (3.44). See Ref.[25] for a related 
problem in the context of a relativistic particle. The canonical Liouville measure is not 
reparametrization invariant in general[26], and it needs to be proved. This problem is 
analogous to the Polyakov-type path integral in string theory[27] - [29]. We here show 
that the naive canonical Liouville measure is gauge invariant up to a renormalization of 
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the cosmological term ( or energy eigenvalue E). This renormalization is however universal 
for any algebraic gauge fixing of the form 



h{r) = f(x(r)) 



(5.3) 



for a gauge fixing function f(x(r)). The analysis presented below, which is formal but is 
known to work in string theory where no simple discretization is known, may be useful 
to understand certain formal aspects such as gauge invariance in the point-particle path 
integral. 

We here analyze the BRST invariant path integral for (5.1) by using the Faddeev- 
Popov procedure, namely, we replace the naive measure in (5.1) by 



/ 



V(Vhx)vVhV(h 3/2 c)VcVB 

xexp{i [ Lhdr + i [ [B(Vh — J~f) — i-cVhd T c — iccd T (Vh — \ff)]dr} 
Jo Jo 2 v 

= J d/j,exp[i J C e ffdr] (5.4) 

The BRST transformation is defined as a translation in the Grassmann parameter , 6 — > 
6 + A, in the superfield notation(note that 6 2 = X 2 = 6\ + \6 = 6c(r) + c(t)6 = 0) 



x{t,6 



y/f(x(r,9) 



c(t,9 
c(t,9 



= x(t) + i9c(T)d T x(T) 

= ^)+i9[c(T)d T + ^(d T c(r))}^) 

= Vhx(r) + i9[c(r)d T + \(d T c(r))]Vhx{T) 

= ] /f(x(r))+i9c(T)d T] /f(x(r)) 

= c(r) + i9c(T)d T c(r) 

= c(t) + 9B(t) 



(5.5) 



For example, the BRST transformation is given by 

5^h[r) = iX[c(T)d T + ^(d T c(r))}^j 

5{Vhx{r)) = i\[c(T)d T + ^(d T c(r))]Vhx(T) 

5(h 3/2 dc(r)) = iX[c(r)d T + l(d T c(r))](h 3/2 dc(T)) 

z 



(5.6) 



In the last relation for the ghost variable c(r), we consider the differential dc{r): one 
may write h 3 ^ 2 dc = d(h 3 ^ 2 c) for a fixed metric h{r), which is the case required to study 
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the BRST invariance of the path integral measure in (5.4). Note that all the variables in 
(5.6) have the same BRST transformation law which is anomaly free [26]; the combination 
VcDB is also manifestly BRST invariant. These properties in turn ensure the BRST 
invariance of the path integral measure in (5.4). One can also confirm that the action in 
(5.4) is also BRST invariant. 

By using the BRST invariance of the action and the measure in (5.4), one can prove 
that the path integral (5.4) is independent of the choice of f(x(r)) as follows: Under an 
infinitesimal change of f(x(r)), the path integral (5.4) changes as 

-i jT dr{[B{T)Sy[J{T) -ic(T)c(T)d T 6yfm]) 

= -iJ dn[B(T)6y/f(ri - ic{r)c{r)d T 5^J{r)}e^o c ^^ (5.7) 

where {B5y/J — iccd T 5y/J) denotes the averaging in the path integral (5.4). Next we note 
the BRST identity 

(c(r)^/(x(r))> 

= J d^{c(r)6^f(x(r))}e l Io^ 

= J d^{c(r)8^f(x(r)) + X[B(r)5^j - ic{T)c{r)d T 5^r)\}e^o ^ 

= (c(T)Sy/f(S(r))) + X(B(r)5^f{r) - ic{r)c{r)d T 5 ^Jij)) (5.8) 

where the primed variables stand for the BRST transformed variables such as \J h{r)' = 
Jh(T)+i\[c{T)d T + \{d T c(T))\y[Hr), 5y/f(x(r))' = 5y/f(x(r))+iXc(T)d T 5y/f(x(r)) and 
c(r) = c(t)+\B(t). The first equality in (5.8) holds since the path integral is independent 
of the naming of path integration variables, and the second equality holds because of 
the BRST invariance of the measure d/j,' = dfi and the effective action C e f/ = Ceff- 
(5.8) shows that (B5-/J — iccd r b\ff) = 0, and thus (5.7) vanishes. Namely, the path 
integral (5.4) is independent of the choice of f(x(r)), provided that one specifies the 
gauge invariant initial and final states. 
The path integral (5.4) is rewritten as 

J VpD£VVhVcDcDBexp{i £ dr[(pS - Hh) + B(Vh- ^) -i^cVhd T (-^c)}} (5.9) 
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with 

H = ^-f + V(r)-E 
2m 

p = Vhp 
x = Vhx 

c = h 3/2 c (5.10) 

In fact, after the path integral over p, one recovers (5.4). We also set y/h — v7 — after 
a partial integration in the ghost sector in (5.9). Note that the path integral 

J Vpexp{i J [--^f + px- l^f\dr} = J Vpexp{i J ^ffidr] ( 5 - n ) 

is a constant independent of the metric h(r): This is another way to see why the variables 
with weight 1/2 such as p and x are chosen for the world scalar quantities p and x as 
reparametrization invariant path integral variables. 

In the above path integral, it is important to recognize that the singular (time- 
derivative) terms in the Lagrangian (5.9) are written as 

/r T 1 1-1 1 

VfiDxVcDcexp{ drHp-^dJ-^x) + ^cVhdJ— ^c)}} (5.12) 
Jo \Jh \Jh 2 n A i z 

Those singular terms have a Weyl invariant structure; namely, the /i-dependence can be 
completely removed by a suitable scale transformation of p, x, c and c such asf-> \/hx 
and c — > (l/y/h)c. In this process of scaling, one obtains a Jacobian (or anomaly), which 
can be integrated to a Wess-Zumino term. At the same time, the path integral variables 
p and x are reduced to the naive ones. In the present case, one can confirm that the 
Jacobian (Weyl anomaly) has a form 

M J hd(3(T)dT (5.13) 

for an infinitesimal scale transformation 

c(r) -> e- <5/3(r) c(r) 
g( r ) _ e M «c(r) 
~p(r) -> eWflr) 

~x(t) -> e^ (T) f(r) (5.14) 
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This evaluation of the Jacobian is performed for p and x, for example, by[28][30] 
Jim / ^ e ~^r J_^uJ_g J_y M 2y £ ikr 

r rlk 1 11 1 

- ^ M Ir^7K WM + *> ve 1[ 7e ( ^ /m + ik) 7n ]) 

r rlk 1 

= J fc «/ s ««p{-* , /* , }(i + o( ]Si )) 

= -^=Mh for M^oo (5.15) 
2V 71 " 

which gives a term of a general structure as in (5.13). The anomaly calculation is specified 
by the basic operators appearing in (5.12) 

Th dr Th ' or (5 - 16) 

both of which give the same form of anomaly proportional to h, as in (5.15). See also 
Ref . [28] . The overall sign of the Jacobian is specified depending on the statistics of each 
variable, i.e., a Grassmann variable gives an extra minus sign. Since only the most singular 
term survives in (5.15), one can confirm that the knowledge of the most singular terms 
in the Lagrangian (5.9), i.e., the time- derivative terms in (5.12), is sufficient to calculate 
anomaly. If one denotes the integrated anomaly ( Wess-Zumino term) by T(h), we have 
from (5.13) 

r 8T r 

T(h) - r(/wT 25/3 ) = J 25(3(T)h(r)-^j-^dr = M J h(T)5(3(r)dr (5.17) 

since h is transformed to he~ 2S/3 by (5.14), and we obtain 

T(h) = M J h(r)dT (5.18) 

with a suitable (infinite) number M. This T(h) is added to the action in (5.9) and it has 
a form of the cosmological (or energy) term in (5.9). 

This calculation of the anomaly is quite general: In two-dimensions, the Weyl anomaly 
has a structure M 2 ^/g + y/gR, [27] - [29]. In one-dimension, the curvature term containing 
R does not exist, and only the cosmological term arises if one performs a reparametrization 
invariant calculation. This Weyl anomaly renormalizes the bare cosmological term in (5.9), 

E + M = E r (5.19) 
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with E r a renormalized energy parameter. But this renormalization is independent of the 
choice of the gauge fixing function f(x(r)) in (5.4). 

We cannot assign a physical significance to the precise value of the renormalization in 
(5.19), since the linear divergence in (5.15) is regularization dependent. Only the Weyl 
invariant structure in (5.12) for the choice of BRST invariant measure and the extraction 
of the metric dependence from the singular terms in (5.12) as a Wess-Zumino term have 
a well-defined physical meaning. 

After this procedure of eliminating the /i-dependence from the time- derivative terms 
and then the integration over B and h, one obtains a path integral 

J VpDxVcDcexp{i J dr[px — f(x(r))H — i-cd T c]} (5.20) 

One may ignore the decoupled ghost sector, which is first order in d T in the present case 
and does not contribute to the Hamiltonian. One thus finally arrives at the phase space 
path integral 

J VpDxexp{i £' dr[p'x - f(x(r))H]} (5.21) 

where the overall normalization factor, which is independent of f(x) but can depend on 
Tf, is fixed to be the same as for f(x) = 1 in (5.21). 

When one cuts the interval [0, 77] into meshes, the path integral measure in (5.21) is 
defined by 

dV f (p) = f(S)d 3 p (5.22) 

for each mesh point, and one performs integral over the momentum first. Since one has 
one extra momentum integration relative to the coordinate one for fixed end points x a (0) 
and Xb{rf), the external states are also specified for the volume element dVf(x). 

A way to provide gauge invariant external states is to consider a trace by integrating 
over x a (0) = x^rf) in (5.21), and let E r — > E n . In this limit the path integral is reduced 
to 

J2 J dV f {x a )<t> n {E n - XaW n (E n ; x a ) = n 2 (5.23) 

where we used the notation in (3.57), and the summation is over the n 2 degenerate states. 
The states (j) n (E n ; x) form a complete physical set for H T (E r = E n ). 
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A less trivial way to incorporate physical external states is to consider the Green's 
function 

G(E r ; x b , x a ) = i J dr f J VpDxexp{i J dr[px - f(x(r))H}} (5.24) 
and look at the pole position at E r ~ E n . In this case, one has (see also eq(3.53) ) 

f dV f (x b )dV f (x a )(pl(E n ]x b )G(E r ]x b ,x a )(p n (E n ]x a ) = ] 

= (E r - E n ) (- d -±0Sf 5 ' 25) 

If one recalls the definition of t/j n (E n ;x) in (3.55), one concludes(see eq.(3.61)) 

G(E r ;x b ,x a ) =Y,^n{E n ;x b ) l _ ip* n {E n ; x a ) (5.26) 

in the sense of an unsubtracted dispersion relation. The pole positions and the residue 
functions in this relation are all gauge independent. 

In the form of the quantity (n\^-\n) in (5.25), one needs to adjust the "wave function 
renormalization factor" \J — dX ^ 3. f or eac h external state to obtain a gauge indepen- 
dent physical quantity, but in the form G(E r ; x b , x a ) in (5.24) one need not supply the 
renormalization factor. 

An interesting implication of the path integral approach(5.4) is that one may under- 
stand the r-integral in (5.24) as an integral over the deformation parameter ( or moduli) 
of an analogue of the world-line. In the present path integral, one sums over all the pos- 
sible deformation of the "world-line" as well as the paths in space to obtain a physically 
meaningful geometrical quantity. Another interpretation of r may be to regard it as an 
analogue of a proper time[31]. 



6 Conclusion 

An attempt to solve the Green's function for the hydrogen atom exactly in the path 
integral [4] opened a new avenue for the path integral treatment of a general separable 
Hamiltonian of Liouville-type. This new view point, which has been shown to be based 
on the Jacobi's principle of least action , provides a more flexible framework of path 
integral to deal with a wider class of problems of physical interest. The Jacobi's principle 
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of least action , besides being reparametrization invariant, gives an attractive geometrical 
picture of particle orbits in a curved space deformed by the potential. On the other hand, 
the fundamental space-time picture of the conventional Feynman path integral, which is 
associated with the Hamilton's principle of stationary action, is lost. ( A 4-dimensional 
picture is however recovered by a generalization of the Jacobi's principle for a relativistic 
particle, S = —m J dr^J (dx^/dr) 2 , and an arbitrary parameter r is identified with a 
proper time in quantum theory [31]). 

In the present paper, we discussed some of the basic issues related to this new approach 
to the path integral from a view point of general gauge theory. We have shown the gauge 
independence of the Green's function. A BRST analysis of the problem from a view 
point of one-dimensional quantum gravity is presented. We also commented on a possible 
application of this scheme to the Stark effect and a resummation of perturbation series. 

I thank C. Bernido for calling the path integral of the hydrogen atom to my attention. 

[Note added] 

After submitting the present paper for publication, the papers quoted in Ref. [32] came to 
my attention. These papers deal with some related matters of the path integral treatment 
of the hydrogen atom. 
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